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SUMMARY 


This paper is concerned with quantifying the effect of sensor and 
actuator errors on the ability of thermal control actuators to correct 
static shape distortion of orbiting spacecraft such as large space antennas. 
Expressions for the statistical properties of the control effectiveness are 
derived based on the assumption of normally distributed errors. 

An example of a 55m space-truss parabolic antenna distorted by orbital 
heating is used to demonstrate the effect of sensor and actuator error. It 
is found that in the case of even moderately large errors, increasing the 
number of actuators does not always enhance shape control. 
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T 

L the Cholesky factor of B (B = LL ) 
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r vector defined by Eq. (6) 

rj jth component of r 

r r vector in the absence of sensor errors 

o 

r r vector based on measured distortion 

m 

T vector of temperature differentials 

T actuator temperature differential in the presence of errors 
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T maximum actuator temperature differential 
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T q actuator temperature differential in the absence of errors 
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u^ structural displacement due to unit temperature increase in ith 

actuator 
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r max 


standard deviation of actuator error 


standard deviation of error in measured distortion 


shape distortion vector 
measured shape distortion 


T 1 / 2 
(t B\|>) 


maximum component of (|» 
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INTRODUCTION 


One of the most important requirements in the design of large space 
antennas is that of surface accuracy [1,2]. While studies have shown that 
under certain conditions high surface accuracies may be maintained by 
passive means [3], as a general rule active controls probably will be 
needed. The disturbances which deform the shape of space structures can be 
divided into two types. One type is transient and, hence, eventually leaves 
the structure unchanged. Transient disturbances can be countered by active 
or passive controls which enhance the damping of the structure. The second 
type of disturbance is typified by those due to manufacturing errors and 
aging [4] which may be considered at least quasi-static if not fixed. These 

latter disturbances may be offset by slowly applied, long-acting 

\ 

corrections. Most research to date has concentrated on controlling the 
transient disturbances by using damping actuators [5]. 
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Interest in the problem of controlling quasi-static disturbances has 
recently increased. Much of the work reported on active control of quasi- 
static disturbances is related to active control of optical systems such as 
mirrors (see [6] for a survey of the state of the art as of 1978). 
Generally, force actuators are employed for static shape control (e.g. [ 7 - 
11]). However, Bushnell [7] characterizes some of these (e.g. [12,133) as 
displacement actuators because they are stiff enough to enforce prescribed 
displacements. Another variation of the force actuator in a truss structure 
is one which changes the length of a member by reeling a cable in or out or 
by using a screw mechanism. This approach is used on some antennas (e.g. 
[14]) to correct fabrication errors, albeit on the ground rather than in 
orbit. A recently proposed alternative is the use of applied temperatures 
on the structure. In [15] the basic concept of using heater actuators was 
described and evaluated on beam and antenna examples. In [ 1 6 ] the problem 
of finding near-optimum actuator locations was investigated and two 
methods based on heuristic integer programming were described and evaluated. 
To date no attempt ha3 been made to study the effects of sensor and actuator 
errors on static shape control. The purpose of the present paper is to 
extend the work described in [15] and [16] to a study of the effect of such 
errors, and to obtain the statistical properties of the control 
effectiveness as a function of the magnitude of the errors. The analysis is 
limited to normally distributed errors and to configurations with a large 
number of sensors. An example of a 55m space-truss parabolic antenna is 
used to demonstrate the effects of sensor and actuator errors. 
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REVIEW OF STATIC SHAPE CONTROL METHOD BASED ON HEATER ACTUATORS 


Although discussed in the context of heater actuators, the methods 
discussed in this paper are applicable to both linear force and temperature 
controls. The equations for temperature control from [15] are briefly 
summarized herein. The reader is referred to [15] for a similar derivation 
for force controls. 

The structure is assumed to be in earth orbit and possess rigid body 
degrees of freedom. The structure is defined over some region n and it is 
assumed that its desired shape has been distorted by an amount t|>(Q) where Q 
is a point in fi and i(i is a vector containing displacement components in 
three orthogonal directions. The distortion is corrected by prescribing 
temperatures at n high-thermal-expansion inserts (actuators) placed in the 
structure. The distortion ij> is assumed to be slowly varying so that the 
actuator inputs may be calculated by a quasi-static analysis. 

The residual displacement 6 is the sum of the shape distortion and the 
correction 


n 

6 (Q) = ><Q) + l u (Q)T (1 ) 

i = 1 


where is the temperature differential of the ith actuator relative to the 

temperature at which ty is measured and u^ is the displacement field due to a 

unit value of T. • 
i 


5 



We seek values of T\ which most effectively offset ^ , that is cause 6 

to be close to zero. A common measure of the smallness of 6 is based on the 
rms value 


2 

^rms 


1 

v 

o 


/ 6 • 6dn 

n 


( 2 ) 


where v q is a reference volume and a dot represents a scalar product. The 
necessary condition for a minimum is 


2 

36 n 

a T “ m3 ' = (2/v ) / (<|» + I u T ) * u i dn *= 0' j = 1,2,...,n (3) 

9T j ° fi i=1 1 1 J 


Equations (3)» a system of n linear algebraic equations for the control 
temperatures, can be written as 


AT = r 


(4) 


where the component a 


ij 


of the matrix A is 


ij 


/ u *u dfl 


and the jth component of the right-hand-side. 


is 


(5) 



The ratio of controlled to uncontrolled rms distortion is called the 
distortion ratio g and is given by 

I 6 • 6 dfl 

g 2 - (7) 

J 

0 tp • rpdd 


When the displacement vector is computed from a discretized model of the 
structure having k degrees of freedom, t|>, u and 6 are finite-dimensional 
vectors of order k instead of vector functions. Equation (1) may then be 
written in a matrix form as 

6 = <|> + UT (8) 

where U is a matrix whose columns are u^, i=1,2,...,n. Equation (2) becomes 

6 Lo = (* + UT) T B (xp + UT) (9) 

rms T 


where B is a symmetric positive definite matrix. Similarly the expressions 
2 

for A, r, and g become 


A «= 


T 

U BU 


( 10 ) 


r 


T 

-U &|> 


(ID 
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( 12 ) 


Qj j + UT) B (\l> + UT) 
T 


1 + 


2ijj T BUT + T T U T BUT 
T 

\|» Bij> 


where use has been made of the symmetry of B. 


EFFECT OF SENSOR AND ACTUATOR ERRORS 


General Equations 

We define T as the vector of control temperatures which would be 
o 

produced in the absence of both sensor and actuator errors. That is, T q is 
calculated from Eq. (*J) in terms of the actual distortion tj; by 


AT 

o 


r 

o 


T 

-U Bi|i 


(13) 


Because of sensor error we measure to instead of to and the control 

m r 

algorithm calculates a corresponding temperature vector as 


AT 


m 


r 


m 



(14) 


Additionally, there are actuator errors (including, possibly computation 
errors) so that the control output is T a - The distortion ratio g is then 

given from Eq. (12) as 
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2 


( 15 ) 


2<J) T BUT + T T U T BUT 
a a a 



Using Eqs. (10) and (13) and the identity T - T + (T - T ), Eq. 15 may be 

a o a o 

written as 


2 

g 



+ 


(T - T ) T A (T - T ) 
a o a o 



( 16 ) 


where g is the ideal distortion ratio (obtained when there are no errors) 
o 


2 2i|i T BUT + T T U T BUT r T T 

-OO o 4 o o 

g o ' 1 + T 1 T 


Statistical Properties of Distortion Correction 

We assume that there are enough sensors to obtain an accurate reading 
of in the absence of errors and that they have no bias. Then f - ^ is a 

random vector with a zero mean. Similarly, we assume that there is no 

actuator bias so that T - T is also a random vector with zero mean. Using 

a m 

Eqs. (13) and (14) we obtain 

T - T = G(* - (18) 

mom 

where 
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_i t 

G - -A U B 


( 19 ) 


We further assume that and T a - T m are independent, normally 

distributed, vectors with covariance matrices and C TT , respectively. 

Based on Eq. (18) T - T is also a normally distributed random vector, and 

mo 

its covariance matrix C TT is 

m 


C TT = E[(T - T ) (T - T ) T ] 
TT m o m o 

m 


GC, G 1 


( 20 ) 


Where E denotes the expected value operator. The total control error T g - 
T q may be written as 


T - T 
a o 


(T " T ) + 
a m 


(T m ” T o) 
m o 


( 21 ) 


and is, therefore, the sum of two independent normally distributed zero-mean 

vectors. Hence T - T is also a normally distributed vector with zero mean, 
a o 

The covariance matrix of T - T is defined as C„„ . From Eq. (21 ) 

a o TT 



C TT + C TT 

m 


( 22 ) 


Next we obtain the expected value and standard deviation of g . For this 
purpose we define the vector V 
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( 23 ) 


T T 1/2 

L U(T a - T o )/(t B\|/) 


T 

where L is the Cholesky factor of B (i.e., B = LL ). Equation (16) may now 
be written in terms of V as 


g 2 - g Q 2 + V T V (24) 

2 2 

The expected value of g , denoted E(g ), and the standard deviation 

2 

o(g ) may be given in terms of the covariance matrix C^ v of V (which from 
Eq. (23) is also a normally distributed vector with zero mean) 


C vv = L T UC tt U T L/* T B* (25) 

a 

From Eq. (24) denoting the components of V as 

E(g ) = g_ + I o (v ) (26) 

i=1 

2 

where a ( v ^) i«e. t the variances of the components of V are the diagonal 
terms of 

- «Vu (27) 


n 



The variance of g 2 is 


2 2 2 2 2 
0 (g ) = E[g - E(g )] 


k 2 2 

E[V V - la (V )] 
i = 1 1 


= E[(V T V) 2 ] - [E o 2 (v )]* 
i = 1 


k k 2 2 k 2 2 

= E E E(v v. ) - [E o (v )] (28) 

i=1 j=1 J i = 1 


It is shown in [17] that 


E(v?v? ) 
i J 


^VV^i^VV^j 


2(C VV>Jj 


using equations (27) and (29) in equation (28) 


k k 

o 2 (g 2 ) = 2E E 
i=1 j = 1 




(29) 


(30) 


While exact expressions are available for the average and standard 
2 

deviation of g , only approximate expressions are obtained for the 
corresponding properties of g. Following the method used in [18] the first 
three terms of a Taylor series are used. Specifically for the square root 
function 
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/ 


1/2 1/2 i _ 1 / 2 1 . 3/2 2 

X =X + -^ X (x-x)--k-x (x-x) 

o 2 o o 8 o o 


(BD 


2 2 

substituting g for x and E(g ) for x we obtain 


- 1/2 2 1 2 2 1/2 2 1 2 22 3/2 2 

g = E (g ) + ^ Lg - E(g ) ]/E (g ) ~ Q [g " E(g )] /E (g ) (32) 


Equation (32) is a good approximation if |g - E(g )| is small compared to 
2 

E(g ). Taking the expected value of both sides of Eq. (32) gives 


E(g) 


(g ) ~ 


c (g ) 

T/2 

8E ' (g 2 ) 


and 

o 2 (g) = E[g - E(g) ] 2 = E(g 2 ) - E 2 (g) 


“ q 2 (g 2 ) _ P H ( g 2 ) 

4 (E(g 2 ) 64E 3 (g 2 ) 


(33) 


(34) 
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ANTENNA REFLECTOR EXAMPLE 


An example of a 55m space-truss parabolic antenna reflector from 
[19] shown in Figure 1 is used to demonstrate the effects of actuator and 
sensor errors on the distortion ratio. The antenna is constructed of 
graphite epoxy while the control elements are aluminum. The reflector is 
distorted from its ideal shape by thermal deformation due to orbital 
heating. The temperature history of the lower and upper surfaces of the 
antenna calculated as in [15] is shown in Figure 2. A previous study [16] 
obtained near optimal actuator locations for controlling the thermal 
distortion with varying number of actuators, and these actuator locations 
were used in the present study. In the present study it was assumed that a 
sensor is located at each of the joints of the truss so that in the absence 
of sensor errors all the displacement components at the joints are known. 
The matrix B of Eq. (9) was taken to be a diagonal matrix with unit entries 
corresponding to the displacement components normal to the reflector 
surface. This means that the mission of the control system was to minimize 
the rms value of the normal displacements at the joints. It was assumed 
that sensor errors are uncorrelated and all have the same standard deviation 
. Similarly, it was assumed that acturator errors were uncorrelated and 

have the same standard deviation o^. 

First, the effects of sensor and actuator errors were studied for 
a 12 actuator configuration from [16] (which had g Q = 0.275). The effect of 

actuator errors on the expected value and standard deviation of the 
distortion ratio is shown in Figure 3. In Figure 3, the actuator 
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temperature error is normalized to the maximum component of the temperature 
vector T . The distortion ratio is normalized to g Q . As may be predicted 

on the basis of equations (25) and (28), E(g) and o(g) vary quadratically 
with a.p. This is due to the actuator temperatures being the solution of a 


minimization problem so that small errors in the temperatures do not have 
any first order effect. However, it is seen from Figure 3 that errors with 
o T greater than about ten percent of the maximum actuator temperature 


differential substantially degrade the effectiveness of the control 
procedure. The bands of width a(g)/g D indicate the scatter of the 

distortion ratio, but since g is not a normal random variable, care must be 
exercised in the interpretation of o(g). Figure 4 shows the effect of 
sensor error on the expected values and standard deviation of the distortion 
ratio. In Figure 4, the sensor error is normalized to the maximum component 
of the distortion vector <Ji. Here, substantial degradation in performance 
occurs for errors which are above ten percent of the maximum component of 
the distortion vector. 

Next the effect of the number of actuators was studied. An error 
combination which has a moderate effect on the 12-actuator configuration was 


selected . 


The value of a, was selected to be 10 
'P 


in. (corresponding to 


0.25 percent of maximum component of ip and the value of o^, was selected to 


be 0.01 degrees (corresponding to 6.5 percent of the 12-actuator maximum 
temperature differential). This combination corresponds to E(g)/g Q of 1.11 

for the 12-actuator configuration. The number of actuators was varied from 

12 to 60 with the same values of and a,. The results are summarized in 

T <p 


Table 1 (a) and also in Fig. 5 which shows the one-o scatter band about 
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E(g) . Figure 5 shows that as the number of actuators increases, the effect 
of the error becomes more severe, to the point that the benefit of 
increasing the number of actuators beyond twenty is questionable. 

The same study was repeated for a smaller error combination and the 
results are summarized in Table 1 (b) and Fig. 6. While it is clear that 
the effect of the error becomes more severe as the number of actuators 
increase, the error is small, and the performance still improves with 
increasing numbers of actuators. 


CONCLUDING REMARKS 


An analytical study was performed to predict and assess the effect of 
actuator and sensor errors on the performance of a shape control procedure 
for flexible space structures using applied temperatures. Approximate 
formulas were derived for the expected value and variance of the rms 
distortion ratio (ratio of rms distortions with and without corrections) 
based on the assumption of zero-mean normally distributed random errors in 
measured distortions and actuator output temperatures. 

Studies were carried out for a 55-meter radiometer antenna reflector 
distorted from its ideal parabolic shape by non-uniform orbital heating. 
The first study consisted of varying the sensor and actuator errors for the 
case of 12 actuators and computing the distortion ratio. It was found that 
substantial inaccuracies occurred when the actuator error exceeded about ten 
percent of the maximum actuator temperature differential, or when the sensor 
error exceeded about ten percent of the maximum component of distortion. In 
a second study, sensor and actuator errors were held constant and the number 
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of actuators was varied. Results of this study indicated that for 
sufficiently large errors, the error magnification, due to increasing the 
number of actuators, can offset the benefits of having more actuators, 
leading to the anomalous result that increasing the number of actuators can 
eventually degrade the performance of the control system. 
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TABLE 1. EFFECT OF SENSOR AND ACTUATOR ERRORS ON DISTORTION 

REDUCTION. g = DISTORTION RATIO NEGLECTING BOTH 
o 

ERRORS, o = ONE STANDARD DEVIATION 


(a) = 10 , o T = .01 


Number of 

S o 

E (g) 

a(g) 

E(g) + a(g) 

E(g) - o(g) 

Actuators 






10 

.3124 

• 3358 

.015 

.3508 

.3208 

20 

.1980 

.2645 

.035 

.2995 

.2295 

30 

.1793 

.2683 

.043 

.3113 

.2253 

40 

.1533 

.2722 

.051 

.3232 

.2212 

50 

.1335 

.2790 

.061 

.3400 

.2180 

60 

.1120 

.2799 

.069 

.3489 

.2109 




(b) 

o=10 , 

o T = .005 


Number of 

g o 

E(g) 

o(g) 

E (g) + o(g) 

E(g) - o(g) 

Actuators 






10 

.3124 

.3185 

.0039 

.3224 

.3146 

20 

.1980 

.2171 

.011 

.2281 

.2061 

30 

.1793 

.2061 

.014 

.2201 

.1921 

40 

.1533 

.1913 

.018 

.2093 

.1733 

50 

.1335 

.1826 

.024 

.2066 

.1586 

60 

.1120 

.1720 

.029 

.2010 

.1430 
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Figure 1.- Tetrahedral truss antenna reflector 
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Figure 2.- Temperature history for antenna reflector 
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Figure 4.- Effect of sensor error on distortion ratio 
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